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Abstract. The combined use of the ESO Very Large 
Telescope (VLT) UT1 Science Verification (SV) images 
and of the Hubble Space Telescope (HST) Hubble Deep 
Field South observations allows us to strengthen the iden- 
tification as a candidate elliptical galaxy of the Extremely 
Red Object HDFS 223251-603910 previously identified by 
us on the basis of NICMOS and Cerro Tololo Interameri- 
can Observatory imaging. The photometry presented here 
includes VLT data in U, B, V, R, I, a STIS unfiltered 
image, NICMOS J, H, and K band data, thus combining 
the 16.5 hours of VLT SV exposures with 101 hours of 
HST observing. The object is detected in all images ex- 
cept the VLT U band and is one of the reddest known 
with B-K= 9.7 ± 0.5. We consider a wide range of models 
with different ages, metallicities, star formation histories 
and dust content, and conclude that the observed spec- 
tral energy distribution agrees best with that of an old 
elliptical galaxy at redshift just below 2. Alternative pos- 
sibilities are discussed in light of their likelihood and of 
the perspective of spectroscopic confirmation. 

Key words: Galaxies: elliptical and lenticular, cD- 
Infrared: galaxies-Galaxies: formation-cosmology: 
observations-early universe 
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1. Introduction 

The formation process of elliptical galaxies is still un- 
known. Direct observations of elliptical galaxies show that 
ellipticals exist up to a redshift of one, but the statistics 
are too poor to draw firm conclusions concerning the evo- 
lution of the population as a whole (see, e.g., Kauffmann 
et al., 1996; Im et al. 1996; and Im and Casertano 1998). 
Indirect evidence based on the evolution with redshift of 
the Fundamental Plane of elliptical galaxies (van Dokkum 
et al. 1998), and on the global correlation of metallicity 
with velocity dispersion (Ziegler & Bender et al. 1997), as 
well as direct measurement of Balmer line strengths (see, 
e.g., Bressan et al. 1996), indicate that the stellar popula- 
tions of ellipticals are old, with ages in excess of 10 Gyrs. 
Unfortunately, the uncertainties in the basic cosmological 
parameters do not allow us to firmly convert age estimates 
into a formation redshift. In particular, if the geometry 
of our Universe is dominated by a cosmological constant, 
stars in ellipticals could have formed at a redshift of 1.65 
and still be older than 10 Gyrs at the present time (with 
n = 0.35, Ov = 0.65, H = 65 km s" 1 Mpc^ 1 ). 

The strong correlation between dynamics and stellar 
population properties argues for a tightly defined forma- 
tion process. On the other hand, models in which ellip- 
tical galaxies form as the result of stellar-gaseous merg- 
ers of disk systems at moderate redshifts appear to re- 
produce successfully the observed properties of elliptical 
galaxies (see, e.g., Kauffmann and Chariot 1998 and refer- 
ences therein); in such models, most stars form before the 
galaxy itself is finally assembled. However, the predictive 
power of these models is limited by our modest knowl- 
edge on star formation processes. In contrast with their 
sophisticated treatment of the evolution of gravitational 
clustering, typical galaxy formation and evolution codes 
compute star formation according to simple empirical re- 
lations which may not be representative of realistic star 
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formation processes. In fact, on the basis of entirely dif- 
ferent theoretical arguments including a detailed descrip- 
tion of ISM properties but no cosmology nor dynamical 
evolution, one could instead argue that spheroidal sys- 
tems are preferentially formed at high redshift (Spaans & 
Carollo 1997, Carollo et al. 1997a). Thus, it appears that 
the present stalemate between different theories and low- 
to-moderate redshift observations awaits an observational 
answer. 

Important clues may indeed be provided by the discov- 
ery of a population of Extremely Red Objects (EROs): 
relatively bright in the NIR (K~ 19-20) and very faint 
in the optical ((R— K)> 5 typically Q). As was imme- 
diately suggested, such a red color can be produced by 
an "old" stellar population at high redshift (e. g., LBDS 
53W091 at z=1.55, Bpinrad et al. 1997) or by a younger 
stellar population enshrouded by dust lanes (e. g., HR10 
at z=1.44, praham fc Dey 1996 , with one detected emis- 
sion line). These objects have generally subarcsecond sizes 
and therefore very little morphological information can be 
gathered with ground based telescopes. Nevertheless, in 
the case of HR10, the morphology appeared to be irreg- 
ular, thus strengthening the case for the object being a 
late-type or irregular galaxy. Furthermore, the presence 
of a young stellar population was suggested also by the 
relatively blue B-I color of HR10 (B-I« 1.7). 

Besides the two prototypical (and particularly bright) 
cases of LBDS 53W091 and HR10 where it was possible to 
determine the redshift using the Keck 10-m telescope, the 
"dust-age" degeneracy in the broad band colors is gener- 
ally coupled with the redshift and metallicity uncertain- 
ties. In other words, the extremely red R-K observed in 
the EROs can be explained by different combinations of 
star formation history, redshift, dust content and chemi- 
cal composition. The only way to photometrically reduce 
this degeneracy is to obtain multiple broad band colors, 
over the widest possible spectral range, and to measure 
morphological properties. 

In this Letter we present photometric information of 
unprecedented completeness about the extremely red i? 1 / 4 
galaxy identified by Treu et al. (1998a) in the Test Image 
of the Hubble Deep Field South. Due to its particular 
location the galaxy has been observed with the VLT-UT1 
and the Hubble Space Telescope (HST), for a total of 9 
bands covering the spectral bands U through K. The VLT 
was able to detect it in the optical down to B=29 ± 0.5 
and to provide us with the optical SED, while HST gives 
us IR colors and multiband morphological information. 

The data strengthen the early identification as an op- 
timal candidate high redshift elliptical galaxy. If spectro- 
scopically confirmed, this identification would imply that 
some ellipticals can be in place already at redshift 2 and 



1 The definition of EROs is not uniform in the literature. 
Other values are, for example, (R— K)> 4.5 (Elston et al. 1991) 
and (R-K)> 6 (Graham & Dey, 1996). 



that their stars have formed at even higher redshifts. In 
Section 2 and the photometric measurements. In Section 
3 we discuss the identification of the object. In this Letter 
the Hubble constant is assumed to be 100 h km/s/Mpc, 
with h = 0.65 where needed. 



2. Observations 

The HDFS NICMOS field was observed with the Test 
Camera on VLT-UT1 during the VLT Science Verification 
(Renzini, 1998) and with NICMOS on board the Hubble 
Space Telescope during the HDFS (Williams et al. 1999) 
campaign. The field was observed for an additional 9 orbits 
with the Space Telescope Imaging Spectrograph (STIS) 
to obtain a broad, high resolution, visible image, and, as 
one of the flanking fields, with the Wide Field and Plan- 
etary Camera 2 (WFPC2). Unfortunately, we were not 
able to include the WFPC2 images in our analysis as, at 
the time of this writing, the final reduced frame had not 
been released. We refer to Fruchter et al. (1999) for the 
description of the observations and the data reduction. 

In the following discussion we shall indicate the Hub- 
ble Space Telscope filters in the following way: C = STIS 
50CCD (clear); Jn = NICMOS F110W; Hi 60 = NICMOS 
F160W; K 222 = NICMOS F222M. 



2.1. Morphology 

The galaxy is clearly resolved in the NICMOS and STIS 
images, as shown in Figure [l]. To strenghten the visual 
classification as an early type galaxy, we fitted the lu- 
minosity profile with an i? 1 / 4 model and an exponential 
model. We performed the fit independently on each STIS 
and NICMOS image using two different techniques: the 
two-dimensional fit code described in Treu et al. (1998b) 
and the luminosity profile code described in Carollo et 
al. (1997b). With the latter technique we fitted also a 
bulge+disk model. We used in each case two different stars 
as PSFs, in order to check the sensitivity of our results to 
a particular choice of PSF. We confirm the Treu et al. 
(1998a) result that only an i? 1 / 4 profile provides an ac- 
ceptable fit. The high resolution STIS observations con- 
firm the size we determined from the earlier NICMOS ob- 
servations; note that the derived half-light radius is larger 
than both the STIS pixel size and PSF FWHM, making 
our result more robust. The uncertainties listed in Table 
[I] include the error on the sky subtraction and the error 
due to PSF estimates. We cannot exclude the presence of 
a faint, extended disk, but we have verified that our re- 
sults are not significantly affected by the presence of such 
disk components by carrying out a variety of bulge+disk 
fitting experiments (e.g., we find that the "bulge" magni- 
tude does not change by more than ~ 0.1 mag). 
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Fig. 1. HST and VLT images of HDFS 223251-603910 C. 
STIS clear . J. NICMOS Jn H. NICMOS Hi 60 K. NIC- 
MOS K 222 . B. VLT B V. VLT V R. VLT R I. VLT I. 

The images are 3'.' 75 on a side. A faint blue source very 
near to the center of the galaxy is visible in the STIS and 
NICMOS J image. 
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Fig. 2. Luminosity profile of the galaxy (squares) in the 
STIS and NICMOS passbands. The PSF-convolved best- 
fit i? 1 / 4 profile is shown for comparison (solid line). 



2.2. Colors 

The multiband colors of HDFS 223251-603910 are sum- 
marized in Table 0. The STIS and NICMOS magnitudes 
are derived by fitting i? 1 / 4 luminosity profiles as described 
in the previous section but were also checked with aper- 
ture measurements and found to be consistent to within 
~0.1 mags. For the VLT data we give aperture magnitudes 
corrected for PSF losses. We estimate the systematic er- 
ror in the correction to be in the range 0.05 — 0.1 mags. 
The uncertainties listed (<5m) do not include systematic 
calibration errors, which might be as high as 0.1 mags for 



the NICMOS magnitudes and as high as 0.05 mags for the 
VLT magnitudes. We adopt here the NICMOS Pipeline ze- 
ropoints. Particularly the Hi 6 o zeropoint used here is 0.08 
mags higher than the one used in Treu et al. (1998a). The 
VLT zero points in the final coadded images, as released on 
the WEB site http://www.hq.eso.org/paranal/sv/, were 



calculated by matching the photometry in the final stacks 
with that of a single photometric night. The overall con- 
sistency of the zeropoints was checked by comparing the 
VLT data with CTIO data and the main WFPC2 HDFS 
field. The colors listed in Table 1 are corrected for galac- 
tic extinction, using the value E(B-V)=0.026 reported by 
Schlegel et al. (1998). 

3. Identification and conclusions 

Using nine broad band magnitudes derived from VLT and 
HST data, we are now able to constrain the SED of ERO 
HDFS 223251-603910 very well. In order to provide a firm 
identification, we computed the colors of a large grid of 
models, with different redshift, age, metallicity and dust 
content, using the population synhtesis code by Bruzual 
& Chariot (1993, GISSEL96 version). At this stage, to 
avoid further degeneracies, we limited our analysis to sin- 
gle burst stellar population models. Different star forma- 
tion histories with the same age and chemical composi- 
tion, would produce intrinsically bluer colors. Specifically 
we considered models with metallicity Z/Z Q =1,0.4,0.2 and 
dust reddening Ay=0, 0.5,1 (using the extinction law from 
Cardelli et al. 1989). Then we computed the x 2 of the 
synthetic colors with the observed ones. In Figure | the 
contour levels of probability 68, 80 and 95 % are shown. In 
the nine panel the metallicity decreases from left to right, 
while the dust reddening increases from top to bottom. As 
can be seen, the redshift is relatively well constrained, but 
the age of the object is much more uncertain. This is a con- 
sequence of the so-called age-metallicity degeneracy (e.g., 
Worthey, 1994), i.e., that broad band colors cannot distin- 
guish the effects of age and metallicity. Only spectroscopy 
could constrain both the metal content and the age of 
the object. The dotted and dashed lines on the upper left 
corner of the plots are the age of Universe as a function 
of redshift (h=0.65; dotted : O = 1, fl\ = 0, dashed : 
fl = 0.3,17a = 0.7). As the galaxy cannot be older than 
the Universe, only the region below the curves is allowed. 
The best fitting models of our grid in the cases A„ = 0, 0.5 
and Z/Z Q = 1,0.4 are shown in Figure ||, together with the 
photometric points. The x 2 and the model characteristics 
are also shown in the figure. 

As a further check we also compared the colors of 
HDFS 223251-603 910 with the co lors of highly reddended 
starbust galaxies ( Calzetti 1998 , and references therein) 
and we did not find any combination that could reproduce 
the measured Spectral Energy Distribution. 

As already noticed in Treu et al. (1998a), the lumi- 
nosity of this galaxy is consistent with its identification 
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Fig. 3. Contour plots of the x 2 . The colors of a grid of 
models (see text for details) have been computed and com- 
pared with the observed ones. The metallicity of the stellar 
populations is Z/Z Q =1,0.4,0.2 from left to right, the dust 
reddening is Ay=0, 0.5,1 from top to bottom. The con- 
tours correspond to a probability of 68, 80 and 95 % from 
the inside to the outside. The two curves in the upper 
left corner represent the age of the Universe as a function 
of redshift for two different cosmologies (see text). The 
object must be younger than the Universe and therefore 
only the region below the curves is allowed. Its redshift is 
better constrained than its age. 

as an elliptical galaxy at z ~ 1.8, e.g., for fl = 0.35 and 
f^A = 0.65, the galaxy would have, after passive dimming 
down to z=0, My = —21.4 and an effective radius of about 
4 kpc. 

Thus, all the photometric data collected so far are con- 
sistent with the identification of this object as a high red- 
shift elliptical. This renders HDFS 223251-603910 a prime 
target for spectroscopic followup. This is within reach with 
IR spectrograph on 8m class telscopes, such as ISAAC on 
the VLT. If confirmed, its identification will demonstrate 
that ellipticals can already be assembled at redshift almost 
2. Furthermore, by measuring its absorption line strengths 
we will determine in great detail its physical properties, 
such as age, metal content and star formation history, thus 
improving our knowledge of the early stages of galaxy for- 
mation. 
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Fig. 4. The best fitting spectra for the cases with 
Z=1,0.4Z© and Ay = 0,0.5 are shown, togheter with the 
observed photometric measurements and the. The filled 
square is the U band VLT upper limit. The error bars on 
the wavelength axis represent the filter bandwidths. The 
main contribuition to STIS CCD clear filter comes from 
the reddest part of the spectrum. The \ 2 with seven de- 
grees of freedom is shown in the upper left corner. 

Table 1. Photometric properties. For every passband to- 
tal magnitudes are listed with their errors, except for U 
and in which case the 95% CL limits is shown. The errors 
on the magnitudes do not include systematic calibration 
errors (see Section ^^). The effective radius, as measured 
by fitting the R 1 / 4 profile to the HST images is also listed. 
The exposure times (texp), the zero points (zp) and colour 
terms (ct), referred to the ESO and STScI released data, 
are shown for reference. 



Band 


m ±5 m 


r c ±Sr c 


tcxp (s) 


zp 


ct 


U 


< 26.5 




17788 


31.39 


0.239(U-B) 


B 


29±0.5 




10200 


33.26 


0.128(B-V) 


V 


26.7±0.3 




14400 


33.87 


-0.066(B-V) 


R 


25.81±0.15 




7200 


34.13 


-0.061(V-R) 


C 


25.38±0.15 


0'.'26±0"06 


25900 


26.15 




I 


24.65±0.15 




10158 


33.55 


0.023(V-I) 


JllO 


22.01±0.10 


O''17±0"04 


108539 


22.16 




Hieo 


20.47±0.10 


0"16±0'.'04 


128441 


21.52 




K222 


19.33±0.15 


0''20±0'.'05 


103163 


20.00 
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